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Given the bases of a matroid, this paper presents a primal algorithm and a dual algorithm for finding the 
circuits of the matroid. 
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1 . Introduction 

A matroid is a combinatorial structure that possesses important combinatorial properties 
of a wide variety of mathematical structures. Such varied structures as vector spaces, trans- 
versals, certain ployhedral corner points, cycles in a graph, spanning trees, and the source arcs 
used by a network flow are all special cases of matroids. Matroid theory provides a convenient 
way to summarize these scattered results and to extend them simultaneously [1, Ch. 21], * 
[2-5], [8], [11]. 

Let E be any finite set of elements. There are several equivalent ways to define a matroid 
on E [1, p. 476], [2], [9], [10]. The following definition is convenient for the purposes of this 
paper. 

Independence Definition of a Matroid Let g be any set of subsets of E with the fol- 
lowing properties: 

(1) If A' ^A and Aeff , then A'e$ and 

(2) for any subset E f of E, all maximal 2 members of g that are contained in E' have the 
same cardinality. 

The members of family ff are called the independent sets of the matroid. The maximal 
members of Q are called the bases of the matroid. It follows from (2) that all bases have the 
same cardinality. Let {g denote the set of all bases. Set {g equivalently defines a matroid on E. 

For example, let E be the set of edges in a connected, undirected graph. The set of all 
forests of this graph form the independent sets of a matroid on E. The bases of this matroid 
are the spanning trees of the graph. Call this matroid a tree matroid. 

A subset of E that is not independent is called dependent. A circuit is any minimal dependent 
set. The circuits of the tree matroid constitute the set of all simple cycles in the graph. 

Let M be a matroid on set E that is specified by its base set © . Form set $ * by taking 
the complements in E of each member of @. The independent sets corresponding to S* also 
satisfy properties (1) and (2) and thus ®* defines a matroid M* on set E. Matroid M* is 
called the dual of matroid M. 

For example, the dual of a tree matroid is the matroid whose bases are the complements 
of spanning trees. Call this matroid the co-tree matroid. The set of all simple cuts forms the 
circuit set of the co-tree matroid. 

What further relationships exist between the bases and circuits of a matroid? 
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Lemma 1 : Let I he any independent set ofmatroid M. 

Let x e E, x i I and I U {x} he dependent. 

Then, I \J {x} contains exactly one circuit. For the proof, see [1, p. 479]. 

Let D be any subset of E. For matroid M, the rank of D, denoted by r(D), is defined as 
the cardinality of the largest independent set of M that is contained in set D. Thus, 
0<r(Z?)<|D|. 

The closure of set D, denoted by c\(D), is defined as the largest set D f such that BCD' 
and r(D)—r(D f ). It can be shown [2, p. 70] that c\(D) is unique. 

Set D is called closed if D=c\(D). For example in a tree matroid, the rank of a subset Z? 
of edges is the cardinality of the largest forest that can be grown inside set D. The closure of 
set D is set D together with all other edges that have both endpoints in the same component 
of D. Consider the tree matroid of the graph in figure 1. Let D={a, b, /}, then, r(Z7)=3, 
cl(D) — {a, b, e,f, g}. Since Dj£cl(D), seti) is not closed. However, set {a, b, e,f, g) is closed. 

A hyperplane of matroid M is any closed set whose rank equals r(E) — l. Let SK denote 
the set of all hyperplanes of matroid M . Consider the tree matroid for the graph shown in 
figure 1. For this matroid, r(E) equals 4 which is the cardinality of each spanning tree (base). 
The set {a, b,f, e, g} is a hyperplane of this matroid since a forest of at most 3 edges can be 
formed from this set. 

Let (?*, £3*, »£?*, (H*, r*, respectively denote the circuit set, base set, independent sets, 
hyperplane set and rank function of matroid M*, the dual of matroid M. 




Figure 1 



Lemma 2: G = [E-H*:H e cH*} 



Proof. Let Ce C Since C is a minimal dependent set of M, it follows that C—{x] is inde- 
pendent for any x e C. Thus, E—C U {x} contains a base of M*, but E—C does not contain 
a base of M*. 

Thus, E—C is a hyperplane of M*. 
Consequently G 6{E-H*: H* e 9{*}. 

Let H* e cH*. Then E—H* contains a base of M. Suppose x f H*. Then H* U [x] is a base 
of M* t and E—H*—{x) is a base of M. By Lemma 1, E—H* contains a unique circuit of M. 
Since E-H*-{x) is independent in M, E-H* is a circuit. Thus G 5_ {E-H*:E* e(H*}. 

Q.E.D. 

Hence to determine the circuits of M we need only determine the hyperplanes of M*. 

2. Algorithms to Construct the Circuits from the Bases 

An algorithm to construct the base set using only knowledge of the circuit set was developed 
by Hull [6]. This section presents two complementary algorithms that construct the circuit 
set using knowledge of only the base set. These algorithms exploit the properties of circuits 
given in Lemmas 1 and 2. 
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Primal Algorithm 

For each Be© and for each x e E, x e B, form the set BU {x}. ~LetB={x u x 2 , . . . , x r(E) }. 
Fori=l, 2, . . . ,r(E), check if B U {#} — {#*} is a base. If so, color x t orange; otherwise, 
color Xt blue. Element x together with all orange elements form a circuit. 

Proof: By Lemma 1, B U {#}, x 4 B, contains exactly one circuit C. If deleting an element 



Xi, i- 



= 1, 2, 



r(2£), creates an independent set B U 



3 - 



{x^, then &< must be present in 



circuit C. Otherwise, x x is not present in circuit C. Consequently, circuit O must consist of x 
and all orange elements. 

Since a cirucit is a minimal dependent set, the deletion of any member x of a circuit C 
results in an independent set /. Thus, there is some base B, I C B, such that C C B U {x}. 
Hence, all circuits are generated by the algorithm. 

Since E is a finite set, the algorithm must terminate after a finite number of steps. 

Q.E.D. 
Improving the Algorithm 

The essential operation of the algorithm is coloring. Coloring involves checking if the set 
B U {x} — {xt} is a base. To generate one circuit requires r(E) such base list checks, i.e., for 
£1, z 2 , . • . , Xrw If the elements of E are treated as letters and the bases written in alphabetic 
order, each base list check is equivalent to checking if a word of r(E) letters appears or does 
not appear in a dictionary consisting only of words of r(E) letters. Of course, computational 
short cuts can be developed for bases with special structures, such as trees. 

Needless to say, it is probable that some circuit C will be generated more than once. 
Repeated generation of the same circuit can be avoided by checking to sec if each set B U {x} 
under consideration contains a circuit previously generated. If so, there is no need to examine 
B [} {x} further since it will yield a previously discovered circuit. 

However, as the number of discovered circuits increases, this becomes more involved. 
Also, this operation is not an identity check as above but a proper subset check which involves 
more computations. Consequently, it might be computationally advantageous to generate a 
circuit repeatedly rather than to perform many proper subset operations. Ultimately, this 
depends on the magnitude of r(E) and the computational efficiency of the proper subset opera- 
tion. See [7, p. 391]. 

Example: Let's generate the circuits of the tree matroid of the graph in figure 2. The bases 
are listed below: 



Given Bases 


(Spanning 
1. 


Trees) 
abde 


2. 

3. 
4. 


abdj 
abej 
acde 


5. 

6. 

7. 


acdj 
acej 
bcde 


8. 
9. 


bedj 
beef 




Figure 2 
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1. Examination of abde 
Consider abde+c 

a — color a orange since bcde is a base 
b — color b orange since acde is a base 
d — color d blue since abce is not a base 
e — color e blue since abed is not a base 

Resulting circuit abc 
Consider abde+f 

a — color a blue since bdef is not a base 
b — color b blue since adej is not a base 
d — color d orange since abef is a base 
e — color e orange since abdj is a base 

Resulting circuit def. 

AH further sets contain either circuit abc or def and will yield no new circuits. 

From section 1 , we know that the circuits of the tree matroid are the simple cycles of the 
graph. Obviously, the circuits abc and def produced by the algorithm are the only simple cycles 
of this graph. 

Dual Algorithm 

Step 1. Determine @ * from @ by taking complements. 

Step 2. For each B* e ® * and for each x e B*, form the set J5*— {x}. 

For each y e E, y tB*, determine if (B*— {x}) U {y} has rank r*(i?) by consulting (B *. 
If so, ignore y. Otherwise, color y red. After examining all y t J5*, let Hb*,x be the set 
B*—{x) along with all red y. For all B* and all x e B*, generate Hb*,x- 

Step 3. Find the complement of each Hb*,x generated in Step 2. These sets are the circuits 
of M. 

Proof: By Lemma 2, we need determine only the hyperplanes of M*. From Step 1, we can 
determine the bases of M*. Clearly, each base B* of M* less one of its members x has dual rank 
r*(E) — 1 and is contained in exactly one hyperplane of M*. This hyperplane Hb*,x is con- 
structed in Step 2. 

Since E is assumed to be finite, the algorithm terminates in a finite number of steps. 

Q.E.D. 
Improving the Algorithm 

Step 1 and Step 3 require only the operation of taking complements. 

Step 2 requires that we test if a set (B*— {x}) U {y} is a base of M*. As before, this re- 
quires checking to verify if this set is on the list of all bases. 

Clearly, it is possible that there exist bases B* x and B* 2 with Xi e B* 1 and x 2 e B* 2 such that 

cl(B^ 1 -{x 1 })=cl(B^ 2 -{x 2 }). 

In this case, the algorithm would generate this hyperplane twice, which is computationally 
inefficient. To avoid generating this hyperplane more than once, the algorithm could maintain 
a list JC of all the hyperplanes already generated. Before the algorithm starts to generate the 
closure of a set 5*— {x} in Step 2, the algorithm could check if set B*— {x} is a subset of any 
hyperplane H* in list ©/?. If so, then J5*— {x} cz H*, and 

r*(j?*~{aj})=r*(j^— l=r*(il*) 

Since closures are unique, it follows that 

H*=Hb*,x- 
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If B*— {x} is not contained in any member of list ©T, then the algorithm should proceed to 
generate Hb*, % as outlined in Step 2. 

Checking list X involves verifying if the set under consideration is a proper subset of any 
member of this list. If the list is long, this can involve a large number of operations. Moreover, 
the proper subset operation is computationally more difficult than an identity check. Conse- 
quently, it might be better to generate a hyperplane repeatedly. Ultimately, this depends on 
the magnitude of r*(E) and the computational efficiency of the proper subset operation. 

Example: Let's generate the circuits of the tree matroid of the graph in figure 2 using the dual 
algorithm. The bases and their complements are listed below: 



Step 1. 





Given Bases 


(Spanning Trees) 


1. 




abde 


2. 




abdj 


3. 




abef 


4. 




acde 


5. 




acdj 


6. 




acej 


7. 




bcde 


8. 




bcdj 


9. 




bcej 


) 2. 


Examination of cf 




1. 


Consider cj—j=c 





Step I Complements 

cj 
ce 
cd 

i 

bd 

af 
ae 
ad 



a — color red since r(ac) = l 
b — color red since r(bc) = 1 
d — ignore d since r{cd)=2 
e — ignore e since r(ce)=2 

Resulting hyperplane abc 

2. Consider cf—c=f 

a — ignore a since r(af)=2 
b — ignore b since r(bf) =2 
d — color red since r(df) = 1 
e — color red since r(ej) = l 

Resulting hyperplane def 

Examination of ce 

1. Consider ce— e=c c is already contained in abc generated above. 

2. Consider ce— c = e e is already contained in def generated above. 

As seen, the examination of each remaining member of @ * will yield no new hyperplanes. 

Step 3. Take the complements of the sets generated in Step 2. 

Set Circuit 

abc dej 

def abc 

3. Primal Versus Dual Algorithm 

The primal algorithm and the dual algorithm are complementary algorithms. 
To generate a circuit, the primal algorithm adds one element to a primal base and then 
deletes elements from the resulting set of r{E)-\-\ members. Each deletion decision requires 
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checking the list of all primal bases. For each primal base, r{E)(\E\—r(E)) such list checks are 
required. Thus, at worst, | £3 \r(E)(\E\ — r(E)) list checks are required in total. 

To generate a dual hyperplane, the dual algorithm deletes one element from a dual base 
and then adds elements to the resulting set of r*(E) — 1 elements. Each addition decision re- 
quires checking the list of all dual bases. For each dual base, r * (E) ( I E\—r *(E)) list checks are 
required. Thus, at worst, | © *\r*(E) (\E\—r*(E)) list checks are required in total. 

Since 

|S|=|S*| 

and 

r*(E) = \E\-r(E), 

both algorithms require the same number of list checks. Moreover, these lists have the same 
number of members. 

Which algorithm should be used? Barring special computational procedures that exploit 
the special nature of the matroid, it is probably best to use the primal algorithm if r(E)<^r*(E) 
and the dual algorithm if r*(E)<Cr(E). In this way, the length of the words on the base list 
will be minimized. 

What if we elect not to generate repeatedly the same circuit or the same hyperplane? 
Which algorithm is best? If the primal algorithm is selected, then we must check if sets of 
r(E)-\-l elements contain any circuit already generated. If the dual algorithm is selected, we 
must check if sets of r*(E) — l elements are contained in any hyperplane already generated. 

If r(E)<jr*{E), the former is preferable computationally. If r*(E)<^r(E), the latter is 
preferable computationally. This is consistent with the above preferences. 

Hull [6] gives an algorithm for constructing the circuits of M from the circuits of M*. If 
r*(£ T )<r(£'), the primal algorithm applied to M* could generate the circuits of M*, and Hull's 
algorithm could be used to generate the circuits of M from the circuits of M*. If r(E)<^r :¥ {E) 1 
the dual algorithm applied to M* could generate the circuits of M*, and Hull's algorithm could 
be used to generate the circuits of M from the circuits of M*. Unfortunately, Hull's algorithm 
is rather involved, and there seems to be no computational advantage in this circuitous 
approach. 
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